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The tenuous atmospheres and ionospheres of comets and outer planet satellites share 
many common properties and features. Such similarities include a strong interaction with 
their outer radiation, fields and particles environs. For comets the interaction is with the 
magnetized solar wind plasma, whereas for satellites the interaction is with the strongly 
magnetized and corotating planetary magnetospheric plasma. For this reason there are 
many common or analogous physical regimes, and many of the same modeling techniques 
are used to interpret remote sensing and in situ measurements in order to study the 
important underlying physical phenomena responsible for their appearances. We present 
here a review of various modeling approaches which are used to elucidate the basic 
properties and processes shaping the energetics and dynamics of these systems which are 
similar in may respects. 

For comets, this interaction produces the often spectacular 
1 INTRODUCTION. ion tail, which is millions of kilometers long. For planetary 

satellites, the Galileo spacecraft and ground-based and Earth- 
Comets and several planetary satellites possess thin orbit based observing have revealed an array of auroral type 

tenuous atmospheres which interact with their outer emissions, 

magnetized plasma environs and solar radiation to ultimately 
drive the atmospheres into escape and producing ions which 
mass-load the outer plasma. In the case of comets the gravity 

is vanishingly small (~ 10 4 times gravity at the Earth's 2. SOLAR WIND INTERACTION WITH COMET 

surface). In the case of the planetary satellites (e.g., Jupiter's ATMOSPHERES AND IONOSPHERES 

largest four moons) gravity is small, but the plasma impact 

processes are strong enough that only a tenuous atmosphere 2.1. The Cometary Atmosphere 

is maintained (see Chapter 1 11.3) and a large neutral escape a well developed cometary atmosphere extends to 

rate, comparable to a moderate comet, is maintained. For distances several orders of magnitude larger than the size of 

these reasons similar approaches are used for modeling the nucleus. Self-consistent models of cometary atmospheres 

various observed phenomena in order to study the underlying must describe phenomena taking place at several spatial and 

causes. The innermost neutral atmospheres have been treated temporal scales ranging from the initiation of the transonic 

with hydrodynamic approaches, and the outermost dusty gas flow near the nucleus (with typical scalelength of 

atmospheres have been treated with the assumption of 10' 2 -1 O' 1 km) to fast chemical reactions (typical scales are 10°- 

vacuum free-flow. In between these two regions is a 10 2 km) to the formation of the cometary ionosphere (typical 

complicated kinetic region where collision mean free paths scales are lOMo 4 km), to the solar wind interaction of the 

are large compared with relevant spatial scales, but not large cometary atmosphere (typical scales are 1 0 5 - 1 0 7 km)., 

enough to ignore. Interactions with cooling infrared radiation Today little direct observational data are available on the 

also become non-LTE because of decoupling between near sur f a ce gas, dust and plasma environments of cometary 
rotational and kinetic temperatures. nuclei. Our present understanding of this region is based on 

Magnetohydrodynamic (MHD) approaches have been used dust coma images from the Giotto spacecraft Keller [1986] 
to study the large scale interactions ol impacting magnetized anc j dusty-gasdynamic calculations having varying levels of 
plasma with the ionospheres of these bodies. For comets, it sophistication [Crifo et al., 1999a]. Verification is primarily 

is the solar wind with its frozen-in magnetic field which based on comparing model predictions hundreds of nucleus 

impacts the cometary ionosphere. For planetary satellites, it radii from the nucleus with in situ and remote observations, 
is the corotating (or nearly corotatmg) plasma flow frozen The b u ik velocity of the outflowing gas is slightly 

into the planet’s strong dipolar magnetic field which impacts subsonic at the surface (the effusion would be sonic without 

the satellite ionosphere. Mass-loading via photo- and charge an y gas-dust interaction: the momentum loss to the entrained 

impact ionization and charge exchange are important for dust particles results in a slightly subsonic outflow), 

explaining the overall structure of both kinds of interactions. Frequent gas-dust collisions erode and drag away grains from 
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the surface and accelerate them to significant velocities within 
a few nucleus radii. Small dust particles (< 10 jim) can be 
fully accommodated to the gas flow and reach velocities up 
to several hundreds of m/s, while larger grains attain much 
slower velocities (the dust terminal velocity is proportional 
to the inverse of the grain size) [Gombosi et al, 1986]. The 
adiabatically expanding gas converts most of its internal 
energy to bulk motion, while it loses some of its momentum 
and energy to the dust flow. This process is expected to 
result in a very cold (several tens of K) hypersonic gas flow 
in the immediate vicinity of the nucleus. The problem is also 
complicated by the fact that dust particles, especially the 
more volatile CHON particle can serve as sources for gaseous 
species. These particles get their name because they are 
composed of mostly compounds from carbon, hydrogen, 
oxygen and nitrogen, as opposed to the more refractory 
silicate dust composed of silicon and metal oxide 
compounds. CHON particles can thus more easily fragment 
and vaporize, modifying the gas and dust mass, momentum 
and energy densities. 

Most of the species observed in cometary atmospheres are 
chemically unstable radicals, which clearly are photo- 
destruction products of more stable parent molecules. 
Chemical reactions between evaporated parent molecules are 
very slow and of secondary importance, while dissociation 
and ionization by solar UV radiation produce highly reactive 
radicals and ions. 

Most neutral coma modeling deals with the two limiting 
conditions which occur far from and near the nucleus: free- 
molecular and fluid-flow models. Free-molecular models 
describe the distributions of molecules in the outer coma 
where intermolecular collisions are rare. These include radial 
outflow models [Haser 1957; Haser 1966; Wallace and 
Miller, 1958; Keller and Meier, 1976], and non-radial 
models accounting for isotropic ejection of daughter species 
[Combi and Delsemme, 1980a; Festou 1981]. Fluid models 
describe the bulk properties of various fluids (gases, dust, 
and plasma) that are collisionally coupled in the inner coma. 
Examples of these are gas-phase chemistry models [Giguere 
and Huebner, 1978; Mitchell et al., 1981; Cochran 1985; 
Allen et al., 1987], combined chemical-dynamic models 
[Huebner and Keady, 1983], dusty-gasdynamic models 
[Marconi and Mendis, 1983, Gombosi et al. 1986], and 
multidimensional dusty-gasdynamic models [Kitamura 
1 986; Komle and Ip, 1987 ; Korosmezey and Gombosi, 1990; 
Crifo 1995; Crifo et al, 1999a; Crifo et al. , 1999b; Combi 
et al, 1999]. 

Monte Carlo methods have bridged the gap between the 
fluid and free-flow regimes in the important transition region 
where many observations occur. After the first Monte Carlo 
coma model [Combi and Delsemme, 1980a], collisions were 
introduced [Kitamura et al, 1985] to model the spatial 
distribution of the H coma and the partial collisional 
thermalization of H atoms. Since then this approach has 
been significantly generalized and extended, making a self- 
consistent picture of the inner and outer coma [Combi 1987; 
Bockelee-Morvan and Crovisier, 1987; Combi and Smyth, 
1988a and 1988b; Ip 1989], added dusty-gas hydrodynamics 
to the inner coma [Combi 1989], a treatment for heavy 
species [Combi et al, 1993], and Ly-a radiative transfer 
[Combi and Feldman, 1992]. These models have been used 


to analyze various observations of water dissociation products 
and trace species in P/Halley including spatial profiles, 
column abundances (production rates), and Doppler line 
profiles. These hybrid models have even been extended to 
treat the large collision region of comet Hale-Bopp (1995 
01). Figure 1 shows outflow velocities and kineitc 
temperatures at a function of distance from the nucleus for a 
number of heliocentric distances as calculated using a hybrid 
ID-spherical hydrodynamics/Monte Carlo calculation [Combi 
et al, 1999]. These predicted outflow velocities and kinetic 
temperatures were shown to provide a reasonable explanation 
for radio observations of the heliocentric distance dependence 
of both outflow velocity and kinetic temperature in the inner 



Distance from the Center of the Nucleus (km) 

Figure 1. Radial Outflow Velocity and Kinetic Temperature in 
the Coma of Hale-Bopp. Shown are the results of a one- 
dimensional spherical hybrid kinetic/dusty-gas 
hydrodynamic calculations from Figure 1 in the paper by 
[Combi et al, 1999]. Above is the variation of the radial 
outflow velocity with distance from the nucleus for different 
heliocentric distances of Comet Hale-Bopp. Below are similar 
curves for the gas kinetic temperature. 
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coma [Biver et al, 1999] and outflow velocity in the outer 
coma [Colom et al, 1999]. 

Fully kinetic iterative test particle (ITP) models devel- 
oped for describing atomic species in planetary exospheres 
have been applied to comets [Hodges 1990]. Xie and 
Mumma [1996a and 1996b] have presented improved ID 
spherical and a 2D axisymmetric versions. In an ITP model 
the trajectories and collisions of molecules of each species are 
computed through a background gas of the modeled species 
themselves, as given by their phase space velocity 
distribution functions. It arrives at steady-state solutions for 
the species distribution functions by iteratively updating the 
background distributions with those accumulated from the 
test particle trajectories. 

In recent years these test particle models have been 
superceded by Direct Simulation Monte Carlo calculations 
(DSMC) orignally pioneered by Bird for high Knudsen 
number gas dynamics flow problems [Bird 1976 and 1994]. 
DSMC is a direct simulation of the generalized rarefied 
problem which is most fundamentally described by a 
solution to the collisional Boltzmann equation. The 
combination of velocity moments of the Boltzmann equation 
with closure arguments and the ideal gas law can be 
manipulated to yield the more familiar (and tractible) 
equations of hydrodynamics: both the Euler and Navier- 
Stokes versions. DSMC-type methods have been applied to 
ionospheric and magnetospheric plasma simulations [Wilson 
et al., 1992; Miller and Combi, 1994, Miller et al, 1995] 
and earth’s corona [Shematovich et al. , 1994]. Pospieszalska 
and Johnson, 1992 and 1995] have presented ID spherical 
particle simulations for Io's S0 2 atmosphere using a method 
which is philosophically between the ITP and DSMC 
methods, mixing elements of both. Most recent DSMC 
calculations for lo’s tenuous atmosphere have been published 
by Marconi et al. [1996] and Austin and Goldstein [2000], 
as we will discuss in a later section. However, of particular 
relevance here is that DSMC has been applied to study high 
Knudsen number flows in neutral cometary atmospheres 
[Combi 1996; Harris et al, 1997; Skorov and Rickman 

1998] . 

Unlike the ion chemistry in the coma which is sparked by 
photochemistry (but where the detailed composition is 
controlled by fast ion-neutral reactions), the neutral chemistry 
is controlled in large measure by photochemistry. The 
composition of a "typical" comet is about 80% water, about 
5 % each of CH 3 OH, and other organics, 3% C0 2 , 0.5% NH 3 
and <1% HCN, CH 3 CN, C 2 H 2 , and C 2 H 6 [cf. Crovisier 

1999] . In addition, there is significant CO and H 2 CO in 
comets (up to 10-20 % total CO in some comets), but the 
fraction that comes directly as a gas from the nucleus, e.g., 
frozen H 2 CO, as opposed to being vaporized from an 
extended source of grains is uncertain. The balance is also 
uncertain between parent molecule gases and vaporizing 
grains as sources for the carbon containing radicals, C 2 , C 3 
and CN. A fraction of the source of CN and C 2 may be 
associated with grains but clearly some of the CN and C 2 
comes from photodissociation of HCN and C 2 H 2 , 
respectively, which are both among the parent gases emitted 
from the nucleus [cf. Brooke et al, 1996; Combi and Fink, 
1997; Rauer et al, 1997]. 


2.2. The Cometary Ionosphere 

The main photoionization product in the cometary coma 
is H 2 0 + . These water ions transfer a proton to other neutral 
molecules, which have a higher proton affinity than the OH 
radical (the prime example is H 2 0 + + H 2 0 -> H 3 0 + + OH). 
The proton affinity of the parent molecules determines the 
reaction rate, therefore protonated ions of molecules with 
high proton affinity are more abundant (relative to other ions) 
than their parent molecules are (relative to other neutrals) 
[Haberli et al, 1997a]. The important species and reactions 
for ion chemistry in a water-dominated coma are summarized 
in tables by Haberli et al [1997a], 

Neutral atoms and molecules of cometary origin move 
along ballistic trajectories and become ionized with a 
characteristic ionization lifetime of 10 5 - 10 7 s. Freshly bom 
ions are accelerated by the motional electric field of the high- 
speed solar wind flow. The ion trajectory is cycloidal, 
resulting from the superposition of gyration and E x B drift. 
The corresponding velocity-space distribution is a ring-beam 
distribution, which has large velocity space gradients and it 
is unstable to the generation of low frequency transverse 
waves. The combination of ambient and self-generated 
magnetic field turbulence pitch-angle scatters the newborn 
ions from the pickup ring. As a result of this process the 
pitch angles of the pickup-ring particles are scattered on the 
spherical velocity space shell around the local solar wind 
velocity [Neugebauer et al, 1991a; Motchmann et al, 1993; 
Coates et al, 1993]. 

Conservation of momentum and energy require that the 
solar wind be decelerated as newly bom charged particles are 
"picked up" by the plasma flow. Continuous deceleration of 
the solar wind flow by mass loading is possible only up to a 
certain point at which the mean molecular weight of the 
plasma particles reaches a critical value. At this point a weak 
shock forms and impulsively decelerates the flow to subsonic 
velocities. The cometosheath is located between the cometary 
shock and the magnetic field free region in the innermost 
coma. The plasma population in the cometosheath is a 
changing mixture of ambient solar wind and particles picked 
up upstream and downstream of the shock. Near the nucleus 
ion-neutral chemistry and recombination starts to become 
more and more important. In the inner coma, an "ionopause" 
or "diamagnetic cavity boundary," separates the solar wind 
controlled magnetized plasma from the magnetic field free 
cometary ionosphere, the diamagnetic cavity. Deep inside the 
diamagnetic cavity the cometary plasma and the neutral gas 
are very strongly coupled by ion-neutral collisions, and they 
move radially outward with the same expansion velocity. 
See Figure 2. 


3. SATELLITE ATMOSPHERES AND IONOSPHERES 
WITHIN PLANETARY MAGNETOSPHERES 

The nature of the atmospheres and ionospheres of Jupiter's 
natural satellites Io and Europa have been treated in chapters 
III. 2 and III. 3. Their interactions with their surrounding 
radiation, and particles and fields environments are very 
active fields of study. Various frameworks, depending on 
different regime-dependent approaches have been adopted in 
recent years, with the hope of understanding the basic global 
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Figure 2. Comet Ionosphere Interaction with the Solar 
Wind. Shown is a plot of plasma flow streamlines (white) 
superimposed on a color plot of Mach number. The 
boundaries of the contact surface and diamagnetic cavity 
are indicated. 

structure of the atmospheres, and their interactions with solar 
radiation and with the Jovian plasma torus environment. 


3.1. Europa's Interaction with Jupiter's Magnetosphere 

Recent theoretical work on the atmosphere/ionosphere of 
Europa has been done using diffusive photochemical and 
energy-balance models. Being constrained by the 
observations of Hall et al. [1995] and Brown [1996], Ip 
[1996] have made ID-spherical calculations that predict a 
lower atmosphere with a moderate exponential scale height 
and an extended corona. Saur et al. [1998] performed a 
study of 3D aspects of the interaction of Jupiter’s plasma 
torus with Europa's atmopshere. They adopted multiple ID 
hydrostatic-type diffusive profiles for the atmosphere and 
account for chemistry and energy balance and the impacting 
plasma torus. They solved for the effects of currents and 
electric fields assuming a constant and uniform magnetic 
field. The recent flybys of Europa by the Galileo orbiter 
indicate the presence of an ionosphere [Kliore et al., 1997], 
roughly consistent with expectations from the O detection, 
and a disturbed local magnetic field [Kivelson et al, 1997; 
Neubauer 1998] which could be accounted for by an intrinsic 
field or an induced field. 

Figures 3 and 4 show the main results of the simulations 
of Galileo’s E4 Europa flyby by Kabin et al. [1999b]. It was 
found that the best agreement with observations of Kivelson 
et al. [1997] and Gurnett et al. [1998] can be obtained by 


assuming a total mass loading rate of 7 x 10 25 0 2 + ions per 
second [Kliore et al,. 1997]. The best agreement with 
observations was obtained when mass-loading was 
constrained to the ram side, assuming that the dominant 
ionization mechanism is electron impact [Kliore et al, 

1997] . Mass loading alone was able to account for the 
observed signature in B z , while reproducing the B* and B y 
signatures required the inclusion of an intrinsic magnetic 
dipole of 65 nT along the y axis. This intrinsic magnetic 
field is consistent with the presence of a subsurface ocean 
[Khurana et al. , 1998; Neubauer 1999; Zimmer et al, 2000] 
which would produce an induced intrinsic field in the (x, y) 
plane with B y being the major component [Khurana et al, 

1998] . However, the induced magnetic field deduced by 
[Khurana et al. , 1998] is somewhat larger than this "best 
fit". This work required adopting a 20° rotation of the 
upstream torus flow toward the direction of Jupiter. A recent 
paper by Paranicas et al. [2000] suggests that this might not 
be a property of the actualy bulk flow of the torus ions but 
simply a local effect such as that which results from the 
electric field rotation in electrodynamic simulations [Saur et 
al, 1998]. 

A later two-species MHD calculation for the plasma 
interaction near Europa has most recently been done by Liu et 
al [2000]. In this calculation the impacting plasma torus 
ions (dominated by 0 ++ , S ++ , and S^) represent one species 
and the local pickup ions (dominated by Ch*) are the other. 
The model parameters and form were otherwise similar to 
those of [Kabin et al, 1999b]. One difference is that with 
the choice of two species there was freedom to be able to 
better estimate electron density because of the different mean 
charge states of the ambient torus ions ( = 2) and the local 
pickup ions ( = 1). Furthermore, the results of the Galileo 
Plasma Science instrument [Paterson et al, 1999] appeared 
in the literature in the intervening time. Figures 5 - 7 show 
the results of the comparison of various Galileo particles and 
fields instrument measurements during the E4 flyby with the 
two- species MHD calculations of Liu et al. [2000]. 

The agreement with the measured electron density data 
from the plasma wave instrument is quite good. Two 
density peaks in the wake region correspond to the closest 
approach and the center of the true corotation wake. These 
peaks consist mainly of C> 2 + . The ambient plasma torus ions 
actually show a small minimum density inside the wake, 
consistent with calculations of flow past an obstacle when 
there is no mass-loading. See, for example, the first reference 
Io model by Combi et al. [1998a], who considered a perfectly 
conducting sphere and no mass-loading. An essentially 
empty wake resulted. The ion speed in Figure 6 shows 
approximately the correct qualitative behavior; however, the 
model produces more of a velocity decrease in the wake. 
Similarly the temperature shows the right qualitative 
structure; however, the model considers total plasma 
temperature which effectively has contributions from both ion 
populations and electrons, whereas the observations are only 
of the ions. In this regime of the plasma torus the total 
electron temperature is uncertain, but it is believed that there 
is an important contribution from a population of high 
energy electrons which is non-negligible toward the total 
plasma pressure. For this reason the plot also compares 
twice the measured ion temperature with the MHD plasma 
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Figure 3. Comparison of measured and simulated 
densities, as well as simulated speed and temperature 
profiles along the Galileo E4 flyby trajectory. 
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Figure 4. Comparison of measured and simulated magnetic 
field vectors (solid lines) along the Galileo E4 flyby 
trajectory. 


temperature. Again the structure with three maxima is 
qualitatively reproduced. 

Lastly, the comparison with the Galileo magnetometer 
measurements in Figure 7 is also quite good. The approach 
of Kabin et al. [1999b] and Liu et ai [2000] has been to 
leave the magnitude and orientation of the internal Europa 
magnetic field as a free parameter and use the comparison 
with the data to converge to the best "fit," rather than 
assuming from the beginning that an internal induced field 
from a sub-surface ocean was responsible. The magnitude 
and direction of the magnetic field in both of these studies, is 
in fact consistent with this picture of the induced internal 
field picture suggested by Khurana et al. [1998] and 
Neubauer [1999], providing an independent confirmation. 


3.2. Io’s Interaction with Jupiter's Magnetosphere 

Our understanding of the interaction of lo with Jupiter’s 
corotating plasma goes back to the discovery of Io-related 
decametric radio emission discovered by [Bigg 1 964] and the 
unipole inductor model which was introduced to explain it 
[Goldreich and Linden-Bell, 1969]. The basic premise of the 
unipole inductor picture has recently been questioned by 
Russell and Huddleson, [2000]. There were a number of 
theoretical studies of this interaction during the immediate 
post-Voyager era (sees [Hill 1983] for a post-Voyager 
review). Early theoretical work was often done either in the 
context of a “thin” atmosphere (e.g., see Cloutier [1978]) 
indicative of the surface temperature ( — 1 30 K), or “thick” 
extended neutral atmosphere (e.g. see Goertz [1980]) more 
indicative of volcanic temperatures (-1000 K). Subsequent 
evidence (see the review by Lellouch [1996]) seems to 
indicate a mixed picture of the global atmosphere, which has 
a large extended corona like a thick atmosphere, but appears 
to be dominated by local major injection of hot (high speed) 
gas/dust plumes to high altitudes but only near active 
volcanic vents. Therefore, although the atmosphere is 
probably only locally thick, it still has a large extended 
neutrat corona which might provide a sufficient source of 
impact ionization and photoionization to explain the plasma 
torus. 

Neubauer [1980] presented an analytical model of the 
Alfven standing wave current system which connects current 
through the ionosphere of lo. Southwood [1980] examined 
data from several Voyager instruments and examined the 
possible role of an intrinsic magnetic field for lo as a way to 
retain a robust enough ionosphere, which could provide 
enough conductivity for completing the lo-Jupiter current 
circuit. Several years after Voyager, 3D numerical studies of 
the plasma flow past lo were performed using electrodynamic 
[Wolf-Gladrow et al 1987] and magnetohydrodynamic 
[Linker 1991] approaches. 

The results of the measurements by the particle and field 
instruments on the Galileo Orbiter during the December 1995 
flyby of lo provided new and important information with 
which realistic simulations for the plasma interaction can be 
tested. The JO flyby occurred nearly in the equator plane of 
lo and perpendicular to a flow wake defined by the corotating 
plasma flow past lo. The spacecraft trajectory passed 
approximately 900 km down-stream of lo (in the sense of the 
plasma torus flow). Along that trajectory physical signatures 
of the wake were seen as a broad depression in the magnetic 
field [Kivelson et al, 1996], sharp peaks in the ion [Frank et 
al, 1996] and electron [Gurnett et al, 1996] densities, a 
slowing of the plasma in the core of the wake, a deep ion 
temperature decrease in the center of the wake, and a large 
(factor of 3) ion temperature rise in the flanks of the wake 
[Frank et al ., 1996]. The magnetic field perturbation was 
broader spatially than the density peak, and showed a double- 
reversed structure, whereby the perturbation (as defined by the 
difference from the outer Jovian B-field value) was actually 
weaker right near the close-approach point than it was 
somewhat adjacent to the center of the wake. 

With the goals to explain the general measured features 
from Galileo and to probe important properties of the 
interaction of the plasma torus with lo’s neutral atmosphere 
and ionosphere, there have been recent efforts to improve and 
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electron density along the Galileo E4 flyby. The dashed 
line at the bottom gives the model density of the upstream 
plasma torus ions, and the dash-dot line gives the modeled 
density of the pickup ions. The electron density is 
constructed assuming the average charge state of the 
ambient torus ions is 2, while that for the pickup 02 + ions is 
one. The model electron density (solid line) with that 
inferred from the plasma wave measurements of Gurnett et 
al. [1998] (open circles). 
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Figure 6. Speed and temperature profiles along the E4 
flyby. The solid lines are the plasma speed (upper panel) 
and plasma temperature (lower panel) from the model 
results. The circles are the ion speed (upper panel and ion 
temperature (lower panel) data of Paterson et al. [1999]. 
Because the single-fluid plasma temperature of the MHD 
calculation is really the sum of the ion and electron 
temperature, which may be roughly comparable in this 
region of the torus, we show the values marked by the stars 
that are twice the measured ion temperature. 

extend the pre-Galileo simulations both in terms of the MHD 
[Combi et al., 1998a; Linker et al. , 1998; Kabin et al., 
2001] and the electrodynamic [Sawr et al. 1999] approaches. 
These two approaches are distinguished by the physical 
assumptions which they each do and do not (or in some 
cases, can and cannot) include. 

The single-fluid MHD approach solves a system of mass, 
momentum, energy, and magnetic field conservation 
equations treating the plasma (ions and electrons) as a single 
fluid. Thus it can treat the deformation (draping) of magnetic 
field lines around the obstacle (the near lo region) and 
produce reasonably realistic flow geometry in three 
dimensions. The imposed outside (upstream) B-field can 
have any value and orientation with respect to the flow 
direction. However, MHD cannot, at least yet, include the 
effects of realistic conductivities (Hall and Pederson) or 
charge separation effects which are likely to be important very 
close to Io where the neutral densities are large and electric 
potentials can introduce non-symmetric flow around the 
body. They either include constant artificial conductivity 
[Linker et al., 1998] or assume perfect conductivity [Combi 
et al, 1998a], however comparisons of the two sets of 
published results do not indicate that this choice has any 
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Figure 7. Magnetic field components along the E4 flyby. The 
solid lines are the results of the model. The circles are the data 
from the magnetometer measurements [Kivelson et al. , 1997], 


important consequence. Both MHD approaches address the 
finite conductivity at the inner boundary, either through 
solutions of Laplace's equation for simple conductivity 
models inside of Io in the case of Linker et al. [1998], or by 
fixing the B-field inside Io to be the time-averaged imposed 
Jovian field by Combi et al. [1998a]. The latter was shown 
previously by Linker et al. [1998] to be the major effect of 
actually solving Laplace’s equation inside the inner boundary. 
The similarity of the results of the so-called fixed-boundary 
model of Combi et al. [1998a] and the so-called conducting 
model of Linker et al. [1998] seems to bear out the validity 
of the simpler approach, as well as the overall consistency of 
the two numerical codes despite their differences in numerical 
algorithms and computational grids. 

Both sets of MHD calculations included mass-loading and 
charge exchange varying spatially as the power law 
distribution found by Schneider et al. [1991] for sodium. 
Linker et al. [1998] presented results for a conducting lo with 
no intrinsic magnetic field as well as for a magnetized Io, 
similar to the vacuum superposition field originally 
suggested by Kivelson et al. [1996]. Kabin et al. [2001] 
considered a somewhat more complicated aspherical neutral 
atmosphere based on a lower atmosphere with an exponential 
scale height distribution [Wong and Smyth, 2000] in addition 
to an extended power-law distribution from the corona, which 
resulted in some minor improvements. Both sets of MHD 
results produce a cold dense plasma wake for nominal levels 
of mass-loading and charge exchange. The non-magnetized 
models produce magnetic field perturbations that are similar 
to the Galileo measurements, but none are quite as deep or as 
broad, and none have the reversal of the pertubation (the 
double-peaked structure) in the center of the wake. The 
magnetized models of Linker et al. [1998] produced a broad 
and deep perturbation, but not the self-reversal at the center of 
the wake (the double-peak or bite-out). Figures 8 and 9 show 
a comparison of the magnetized model of Linker et al. [1998] 
and the non-magnetized models of Combi et al. [1998a] with 
the measurements. They do not enable a definitive statement 
to be made regarding an intrinsic field. A recent review of 
the later Io flybys still yields similarly ambiguous results 
[Kivelson et al., 2000] in this regard, although there is some 
hope of a conclusive answer from a successful polar flyby 
during the final Galileo extended mission. 
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None of the standard MHD models have been able to 
reproduce the factor of three temperature rise in the flanks of 
the wake. This temperature rise is roughly coincident 
spatially with the broadened extent of the magnetic field 
perturbation. It was shown by Khurana et al. [1997b] that if 
the temperature increase were simply added in an ad hoc 
manner, it would account for all or most of the broadened 
extent of the magnetic field perturbation. In addition to 
including a more complicated neutral atmosphere distribution 
(as already mentioned), Kabin et al [2001] also included 
deposition of pickup ions with their local corotation plus 
gyration energy in order to explain the broad magnetic field 
disturbance and the large temperatures in the flanks. The 
results generally verity the contention of Khurana et al, 
1997b] but do not address the source of the extra energy. It is 
noteworthy to point out that the simulations of Linker et al. 
[1998] all assume that the ambient (upstream) plasma 
temperature is 200 eV, whereas the calculations of Combi et 
al. [1998a] adopted a plasma temperature 92 eV from the 
measured ion results of Frank et al [1996]. A careful 
comparison of the measured temperatures to the models of 
Linker et al. [1998] shows that they actually have the same 
small 10% increase of temperature in the flanks of the wake 
as the models of Combi et al. [1998a]. This increase is 
explained by the injection of nearly stationary pick-up ions 
into the supersonic plasma flow [Szego et al, 2000]. The 
model temperature profile of Combi et al. [1998a] matches 
the Galileo data far from the wake (92 eV) and in the center 
of the wake, but does not increase to 300 eV in the flanks of 
the wake. The model temperature profile of Linker et al 
[1998], on the other hand, is a factor of two too large away 
from the wake, cuts through the increase in the flanks, and 
then matches the cool observed temperatures in the center of 
the wake. 

The electrodynamic approach Saur et al. [1999] starts from 
the assumptions that the magnetic field is everywhere 
uniform with the local Jupiter dipole value at lo, and that the 
flow is initially perpendicular to this field. The electric 
potential equation is solved in a series of planar cuts 
perpendicular to the magnetic field and through Io. The 
resulting flow is restricted from traveling along the magnetic 
field direction although it can be deflected around the 
obstacle within the plane from which it started. Hall and 
Pederson contributions to the conductivity are calculated 
using a neutral atmosphere having a scale height distribution 
similar to their previous work for Europa [Saur et al, 1998]. 
Their conductivity, however, goes to zero at a small altitude 
above the surface. They also include separate descriptions for 
ions and electrons. From the electnc field described by the 
conductivity and charge separation, they calculated the 
deflected flow and the magnitude of the perturbation to the 
uniform B-field. Their source terms are self-consistent, rather 
than fixed as in the MHD models, however, at their present 
level of their model sophistication the magnetic field 
perturbation could not be fed-back into the plasma trajectory 
calculations self-consistent ly, so that they do not account for 
where the magnetic field is most strongly disturbed. Saur et 
al. [1999] performed many one-dimensional path-integral 
calculations along trajectories defined by the electric potential 
solution to obtain a full three-dimensional picture of ion and 
electron density, velocity, and temperature. 
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Figure 8. Comparison of the Mass-Loading Models with 
Galileo Particle and Field Measurements. The magnetic field 
measurements in (a) from Kivelson et a l [1996] and the 
plasma density (b), pressure (c), and temperature (d) from 
Frank et al [1996 obtained from the plasma ion instrument are 
shown with the Combi et al. [1998a] model values. The solid 
lines give the results for the model with mass-loading and 
reflective boundary conditions. The dashed lines give the 
results for the model with mass-loading and fixed boundary 
conditions. 


Therefore, both the MHD and electrodynamic approaches 
have their strong and weak points. The high conductivity 
assumptions of the MHD calculations are more appropriate 
farther away from Io and the detailed conductivity 
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Magnetized Io 



Time 

Figure 9. Comparison of simulated data traces from MHD 
computations of plasma flow past Io (gray curves) with the 
Galileo data (black curves). Results are from the magnetized 
model of Linker et al. [1998]. Data from the simulations 
were extracted along the Galileo trajectory for the Io flyby. 

(a) Fractional change in the magnitude of the magnetic 
field. Note the “double-valley’ 1 structure in the Galileo data. 

(b) Plasma density, (c) Velocity magnitude, and (d) 
Temperature. 


calculations of the electrodynamic approach are better closer 
to Io's surface. Another important difference is that because 


the ion motion is restricted to the plane perpendicular to the 
average Jovian magnetic field, the flow in the electrodynamic 
approach is essentially incompressible (and inexpandable) and 
therefore cannot collapse in the wake behind lo as 
accomplished in MHD. Therefore, whereas the MHD 
calculations form a dense cold wake similar to the Galileo 
measurements, the approach of Saur et al. [1999] produces an 
empty wake. 

In follow-up work published in his Ph.D. thesis, Saur 
[2000] shows that the combination of their assumed uniform 
B- field with their calculated perturbation field, can produce a 
double-peaked structure in their calculated magnetic field. 
However, this extra bite-out coincides with the location of 
the empty wake, and it is not clear that the extra bite-out 
would be present if their calculation produced the observed 
cold, dense wake as is observed and as the MHD calculations 
reproduce. Unfortunately, what was stated by Linker et al. 
[1998] remains true: no simulation to date has produced 
quantitatively and simultaneously all features of the Galileo 
measurements: the dense, cold plasma wake, the broad extent 
and the central bite-out in the B- field across the wake, and 
the large temperature peaks in the flanks of the wake. Clearly 
simulation science has a long way to go before being able to 
address self-consistently the physically realistic complex 
nature of this coupled system. In the mean time, while 
acknowledging the strengths and weaknesses of both these 
and other approaches, much can still be gained by further 
studies using all of these approaches, which are complemen- 
tary in many ways. 

The electrodynamic approach could be generalized to 
include a realistic nonuniform magnetic field. This, among 
other things, would permit the flow to collapse behind the 
satellite and hopefully produce a realistic dense wake. The 
MHD approach could be generalized to include multiple 
species to account separately for newly deposited pickup ions 
and the background torus ions. Another possible refinement 
would be to include parallel and perpendicular temperatures 
(and pressures). It is quite clear that the newly deposited 
ions, either through mass-loading by new ions or by charge 
exchange, are introduced with a ring beam distribution. Both 
MHD and electrodynamic approaches need to better address 
the range of realistic conductivity around Io. In principle the 
MHD approaches treat the plasma to be perfectly conducting 
[Combi et al., 1998a]. The equations of Linker et al. [1998] 
showed an ohmic heating term, however it was not included 
in most of the simulations. They did include extra viscous 
terms to capture sharp fronts. The electrodynamic approach, 
on the other hand, assumes more realistic Hall and Pederson 
conductivities at low altitudes. This form of the conductivity 
goes to zero at finite altitude, rather than asymptotically 
approaching large conductivities dominated by ion-electron 
collisions at high altitudes [So/ig et al, 2001]. 


3.3. Ideal MHD for Gyration Dominated Plasma at Io 

As discussed above, recent MHD models [Combi et al, 
1998a; Linker et al, 1998] for the interaction of the Jovian 
plasma torus with the tenuous atmosphere of Io reproduced 
many aspects of the December 1995 Galileo Io fly-by 
measurements in the wake, such as the general flow 
geometry, the velocity perturbation, the density peak and the 
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temperature trough [Frank et al., 1996; Gurnett et al., 
1996b]. Only the intrinsic field model of Linker et al. 
[1998] produced a magnetic field signature approaching the 
breadth measured [Kivelson et al, [1996], and neither 
reproduced the measured temperature peaks in the flanks on 
either side of the wake [Frank et al., 1996], which exceeded 
the nearby - 92 eV torus values by a factor of 3. Combi et al 
[1998a] concluded that a process which adds energy in the 
flanks of the wake and/or upstream of the wake was required 
that was not included in the published MHD calculations. 
This was based in part on the suggestion by Khurana et al 
[1997b] that if a temperature increase in the flanks were 
somehow added, a broader magnetic field disturbance would 
naturally result. 

There are some important issues regarding the structure of 
the neutral atmosphere which have been published since the 
previously published MHD models [Combi et al, 1998a; 
Linker et al, 1998], and so a new look is warranted. In 
addition, we have reconsidered one of the basic model 
parameters assumed before, the ratio of the specific heats, y, 
which was chosen to be 5/3. It is interesting to note that 
often in cometary and solar-wind pickup ion simulations, 
investigators have noted that in situations where the pickup 
ions remain in pickup ring distributions, it can be argued that 
there are really only two translational degrees of freedom 
because two are effectively tied together [Flammer et al, 
1991; Ellison et al, 1999]. In this case, the appropriate 
value for y is 2 instead of 5/3. It is quite clear from the 
Galileo magnetometer data that the ion-cyclotron wave 
signatures of pickup ions in the wake of lo are present. This 
is caused by the fact that neither collisions nor wave-particle 
instabilities are efficient enough to isotropize the distribution 
functions of particles during the time-scale of the plasma 
interaction near Io. In fact, the presence of a wake footprint 
for Io seen in Jovian aurorae [Clarke et al, 1998] more than 
half-way around Jupiter is an indication that particles may 
remain in a ring-beam distribution far downstream from Io. 
Although we might describe the background plasma with 
some thermal energy, it is clear that the nature of any 
magnetized plasma interaction is more gyration-like than 
isotropic in nature. Therefore, following the lead of other 
problems we explore here the effect of using a value for y of 

2 instead of 5/3. 

For the first set of new model calculations presented in this 
paper the parameters remain the same as in the previous work 
of Combi et al [1998a], except for taking y to be 2, and its 
effect on the Mach number which now also becomes 2.0 
instead of 2.2. The purpose here is not to present a 
substantially new step, but rather to combine a few features 
from published work to explore possible improvements in an 
up-to-date MHD calculation. For this case we use the same r 

3 5 spatial distribution for mass- loading and charge exchange 
as for both previous MHD models [Combi et al, 1998a; 
Linker et al, 1998], where r is the distance from the center of 
Io. The remaining model parameters are: 

upstream plasma density = 3500 cm 3 

upstream plasma temperature = 92 eV 

upstream mean molecular mass = 22 amu 


upstream magnetic field 

= 1800 nT 

corotation flow speed 

= 56.8 km s’ 

Alfven Mach number 

- 0.4 

Mach number 

= 2.0 

ratio of specific heats (y) 

= 2 

We use the fixed boundary 

conditions at Io which 


correspond to setting the plasma density, velocity and 
magnetic field in ghost cells at the inner boundary. A 
physically-consistent boundary condition either at Io (or at 
Jupiter for that matter) relates to the magnetic field 
perturbation via finite ionospheric conductivity. We must set 
values of plasma velocity and density and magnetic field at 
the boundary to achieve the desired result. As mentioned 
above, the fixing of an internal B-field at the average Jovian 
field level at Io’s location produces the same results as 
assuming a constant finite conductivity inside the inner (Io) 
boundary in the explicit calculations of Linker [1991]. We 
again set the lower Io boundary at an altitude of 1 50 km, the 
plasma density at 10500 cm* 3 , and the magnetic field at the 
time average value from Jupiter’s dipole field at Io (1800 nT). 
In their recent paper Linker et al [1998] performed a similar 
calculation with a somewhat more complicated layered 
conductivity model inside lo. However, the similarity 
between their conducting model and the fixed boundary 
model (not including the different assumed upstream torus 
plasma temperatures) indicates that the conditions at the 
lower boundary between the two MHD groups are in fact 
reasonably consistent with one another. Finally, as before 
the simulation volume for the new calculation is 900 x 600 x 
600 Rj 0} and there are 9 levels of refinement and 92,000 cells 
with sizes ranging from 0.1 R lo near Io to 50 R| C far upstream 
and downstream. 

The change of y to two seems to provide a natural 
mechanism to explain at least some of the previous model- 
measurement differences, namely the temperature peaks on 
either side of the flanks of wake. The previous fixed- 
boundary model produced not only a density peak and 
temperature drop in the center of the wake and the general 
velocity field, but also a deeper magnetic field disturbance 
and a large increase in the temperature in the flanks of the 
wake, as shown by the thin solid lines in Figure 10. Also 
shown are the previous results [Combi et al, 1998a] 
indicated by the dashed lines. The mass-loading and charge- 
exchange rates are the same as the previous result. We make 
no attempt to improved the fit of the density peak by 
lowering the mass-loading but simply show what the effect is 
of changing y. 


3.4. Mass-Loading by a More Realistic Neutral Atmosphere 

We also address here the spatial distribution of the mass 
loading and friction terms by incorporating a more realistic 
neutral atmosphere distribution. Essentially this combines 
the most realistic features of the extended neutral corona 
adopted in the MHD simulations [Combi et al, 1998a, 
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Linker el al, 1998] with the scale height distribution of in 
lower atmosphere used in the electrodynamic simulation of 
[Saur el al ., 1999] and some recent atmosphere models. 

The r' 3 5 spatial distribution for the ion mass-loading and 
charge exchange (friction) terms, which was adopted for the 
MHD simulations from the sodium eclipse measurements 
[Schneider et al., 1991], underestimates the expected rise in 
neutral density near Io. Models for the neutral atmosphere 
consist of a dense surface atmosphere near the surface 
temperature (115-130K), and the corresponding scale height 
of only ~ 12 km [Summers el al., 1996; Wong and Johnson, 
1996, Wong: and Smyth, 2000]. In this region (which MHD 
models do not resolve anyway) there is not substantial 
ionization. At increasing altitudes, of importance for 
understanding the plasma interaction, the scale height 
increases to eventually to values on the order of 100-200 km. 
At still larger distances a partially escaping corona forms 
with a power-law distribution like the sodium distribution as 
modeled by Smyth and Combi [1997]. To model the overall 
plasma interaction, a spatial profile was adopted which is 
consistent with these neutral atmosphere measurements and 
models. In particular, a scale height form for the near-surface 
atmosphere was chosen to be relevant for the altitudes that 
can be resolved by the computation grid and that merge into 
a power-law for higher altitudes. 

In the region from 1.4 to 7 Ri 0 a power-law distribution, 
varying as r 3 5 was found for the observations of the neutral 
sodium corona by Schneider el al. [1991]. However, because 
the bulk of the neutral corona does not consist of sodium^ 
but rather the major SO 2 byproducts, S and O, the earlier r 
distribution may not be appropriate for the bulk constituents. 
The models of Smyth and Combi [1997], which reproduce 
the sodium corona measurements of Schneider et al [1991], 
correspond to the non-thermal kinetic escape mechanism of 
Sieveka and Johnson [1984] with an most probable upward 
flux velocity of 0.5 km s' 1 . Sodium has an ionization 
potential of only 5.1 eV, and in typical plasma torus 
conditions has a lifetime against electron impact ionization 
that is a factor of 2 to 8 shorter than for S and O, which have 
ionization potentials of 10.36 and 13.62 eV, respectively. 
Accounting for the difference in the relative loss of upwardly 
moving Na atoms compared with the average of S and O 
implies that the bulk corona density (as opposed to just the 
trace amounts of neutral sodium) should decrease with 
increasing distance as r 2 5 rather than the steeper r 3 5 for Na. 
Such a distribution is also consistent with some very recent 
Hubble Space Telescope measurements by Wolven et al. 
[2001, in press] of the UV emission by atomic S and O. 

Therefore, we have adopted a new form for the radial 
dependence of the mass- loading and charge-exchange rates to 
be A exp(-z/H) + B r' 25 , where A and B are constants 
adjusted during model-data comparison, z is the altitude, r is 
the distance from the center of Io, and H is an effective scale 
height set to 1 14 km for the region in question. As before, 
the corona is cut-off at 7 Io radii from the center of Io. We 
adjust only one extra free parameter compared with previous 
MHD models, the relative contributions from the lower 
atmosphere and the extended corona, with the transition 
being around 2600 km from the center of Io, and loosely 
corresponding to the atmospheric exobase. Ther inner Io 
boundary has been moved from 1 50 km above Io down to the 


surface, and three more levels of refinement were added 
making the size of the computational cells about 22 km near 
the surface and raising the total number of cells to about 
150,000. 

In Figure 10 the Galileo PLS and MAG measurements are 
compared with the fixed-boundary model of Combi et al 
[1998a], that model but with 7 = 2, and the two-component 
atmosphere model with y = 2. The two-component model 
reproduces the combination of magnetic field, plasma density 
and temperature better than either published MHD modeling 
attempts. The somewhat elevated temperatures in the flanks 
of the wake result from a combination of the use of y = 2 
and the adjustment of the lower and upper atmosphere 
components. As before the velocity flow field does not 
change appreciably between the models, which appears 
essentially the same as shown by Combi et al [1998a]. 

The total fresh ion mass-loading rate and the charge 
exchange rate, which contributes to the momentum and 
energy friction terms, is 1.2 X 10 28 s 1 . The relative 
contributions from the inner atmosphere and outer corona axe 
in the ratio of 15 to 1. The new outer corona contribution is 
smaller than the previous model results. The velocity of ions 
is again similar to previous results, the variation of which 
matched the measurements by Frank et a l [1996]. The 
asymptotic velocity increase to corotation down the wake is 
similar to the values from the radio occultation measurements 
[Hinson el al, 1998] obtained in July 1997. These imply 
that the flow in the wake accelerates from near 30 km s' at 
around 2 Rj c to nearly corotational speeds by only 5 to 7 Ri 0 
down the wake. On the other hand, observations of the 
footprint of the wake of Io by Clarke et al [1998], which is 
weakly present even half-way around the orbit of Io (hundreds 
of Rjo), have been interpreted as indicating a persistence of 
subcorotational flow far down the wake. 

The reversal of the magnetic field disturbance at the center 
of the wake [Kivelson et al, 1996] is still not reproduced; 
the magnetic field disturbance in B z is still not broad 
enough; and the temperature peaks in the flanks are not 
extended spatially enough. However, with the modified 
neutral atmosphere distribution, the new model results do 
provide small temperature peaks and at least an inflection 
point and a flat distribution in B z at the right locations. 
There are a number of possible explanations, which new 
measurements and modeling studies might provide. The first 
would be an intrinsic magnetic field at Io. The 124 flyby of 
Io [Linker et al, 1999a] is inconclusive regarding an intrinsic 
field, and unfortunately the failure of the 125 flyby will not 
yield the anticipated conclusive measurements. It is possible 
that the recently announced flyby in the last extended Galileo 
mission phase may eliminate that as a viable explanation. 
The spatial characteristics of the mass- loading and charge 
exchange source rates may just be complicated. None of the 
MHD calculations include a self-consistent ionization source 
which is based on electron impact ionization of a given 
neutral atmosphere. More realistic simulations could also 
involve latitudinal and/or longitudinal variations, or irregular 
(volcanic) atmospheric distributions. Finally, the inclusion 
of other physical processes in simulations, such as those 
associated with the expected pressure anisotropy, multi- 
species effects, electron beams, currents, other non- 
equilibrium or kinetic effects, or other energy sources [Kabin 
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Figure 10. Comparison of lo MHD models results with 
Galileo measurements. Shown is a comparison of Galileo 
PLS [Frank et al., 1996] data and MAG [Kivelson et al, 
1996] (the open circles) with originally published (dashed 
lines) results [ Combi et al, 1998a], the same spatial model 
with y " 2 (thin solid lines), and a two-component 
atmosphere model with y=2 (thick solid lines). The two 
top panels give the plasma density and temperature, and the 
bottom three give the magnetic field components. 

et al, 2001; Khurana el al., 1997b] might be required to 
substantially improve the agreement between models and 
data. 


3.5. The Interaction of Titan with its Plasma Environs 
Titan possesses a dense local atmosphere [Kuiper, 1944; 
Lindal et al., 1983] which in the presence of an exterior 
magnetized plasma flow necessarily results in a strong and 
complicated interaction. Significant efforts have been 
invested into modeling Titan’s upper atmosphere and 
ionosphere. These are discussed in chapters I.I., 1.2 and 1.3 
of this monograph. Titan is also an interesting because it 
spends a large fraction of its orbit inside the Saturnian 
magnetosphere and a smaller part outside it, in the solar 
wind. This occurs when the solar wind flux is large and the 


Saturnian bow shock is pushed closer to the planet. 
Therefore, the plasma environment around Titan can vary 
considerably. During the Voyager 1 encounter, Titan was 
clearly inside Saturn’s magnetosphere. 

Voyager 1 was the only spacecraft to visit the vicinity of 
Titan, having a closest approach on 12 November 1980 of 
2.7 Titan radii from the center of Titan. These are the only 
magnetic field and plasma data available to the scientific 
community until Cassini has its first flyby of Titan in 2004. 
These measurements show a strong interaction between Titan 
and Saturn’s magnetospheric plasma [Ate et al,. 1982; 
Hartle et al,. 1982]. The measurements and later analysis 
[Neubauer et al, 1984] suggest that Titan’s magnetotail was 
produced by the draping of corotating magnetospheric field 
lines about Titan’s ionosphere. The Satumward aberration of 
the wake was ascribed to a deviation of the corotational flow 
by about 20° from tangential. An asymmetry in the magnetic 
field strength was explained by the fact that the solar 
direction was at a large angle from the incident plasma flow, 
producing an inbound/outbound difference coincident with 
the dayside/nightside hemispheres. Neubauer et al. [1984] 
suggested that the difference is related to two different 
magnetic field-line draping regimes associated with light and 
heavy ions. 

A number of comparisons were made between Titan and 
the classical induced (comet-like) magnetosphere of Venus 
[Kivelson and Russell, 1983; Veregin et al, 1984; Luhmann 
et al, 1991]. While the wake structure of Titan is similar to 
that of Venus, it must be remembered that the flow upstream 
of Titan is sub-Alfvenic in Saturn’s magnetosphere whereas it 
is always super- Alfvenic at Venus. In the last few years, 
encouraged by computational and numerical advances and by 
the interest renewed by the upcoming Cassini missions a 
number of numerical simulation studies have been 
undertaken. 

Modeling Titan’s interaction with the Saturn magneto- 
sphere, or the solar wind, is complicated by the fact that the 
validity of a fluid approach is questionable. Fluid 
descriptions (like ideal MHD) are only rigorously applicable 
in situations where the conductivity is high, the plasma 
parameter (number of charged particles in a Debye sphere) is 
large, the ion gyroradius is small compared to the typical 
length scale of the problem and characteristic speeds (sonic 
speed and Alfven speed) are much less than the speed of 
light. For Titan the only of these criteria that possibly 
becomes questionable is the gyroradius, which for a thermal 
N + ion is 2.25 R-ntan- For the much lighter H + , the 
gyroradius is only 0.16 RTuan- Far from Titan, this is not an 
issue, as tong as computational cells are larger than 
gyroradius. However, close to Titan, both the plasma speed 
and temperature decrease (at - 1.2 Riit an they are about a 
factor of four smaller than upstream values) while the 
magnetic field intensity is about a factor of two larger near 
the front stagnation point. This yields about an order of 
magntitude decrease in the gyroradius near the planet, 
keeping the conditions closer to the realm of ideal MHD. 
Luhmann [1996] further suggested that wave-particle 
interactions and plasma instabilities would also tend to 
decrease gyroradius. 

Cravens et al. [ 1 998] produced very useful three-species 2- 
D calculations of the plasma flow past Titan. The calculation 
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successfully reproduced several main features of Titan’s 
magnetospheric interaction, however, as in any 2-D model, it 
missed a number of important essentially 3-D effects. 
Ledvina and Cravens [1998] presented the first 3-D single- 
fluid calculations, where Titan was approximated as a region 
of high plasma density surrounded by an ion-producing shell. 
They used a nonuniform structured Cartesian grid with the 
small cell size of 0 . 1 R T!la n at the surface of Titan and with 
upstream conditions imposed at 10 R-ntan from the satellite. 
Their simulations produced a narrow wake and Alfven wings, 
which are features detected by Voyager or expected from 
theoretical consideration. This and subsequents simulations 
found ion mass- loading rates on the order of 10 s . 
Intercomparison of variations between models, in this regard, 
is complicated by differences in the models’ inner boundary 
conditions, computational cell resolution, and number and 
classes of species. 

Kabin et al [1999] presented results for a 3-D single-fluid 
calculation making use of their adaptive grid numerical 
scheme which simultaneously included finer resolution near 
Titan as well as a much larger simulation domain. Their 
simulations included a well resolved central body and the 
better numerical method and computational grid was able to 
reveal more detailed structures in the magnetospheric 
interaction. Kabin et al. [2000] introduced a model with an 
internal magnetic field, which they suggested represented the 
induced field created by the ionosphere and which was able to 
better represent asymmetric aspects of the Voyager data. 
They also considered the effects of asymmetric mass-loading. 
Subequently, Kopp et al. [2001] argued that nonsymmetric 
mass-loading similar to that of Kabin et al. [2000] can be 
used to model the effects of the large ion gyroradius in a 
fluid description. The model of Kopp et al. [2001] provides 
an interesting connection between fluid and kinetic 
approaches. Another example of a combination of fluid and 
particle description is a work by Ledvina et al. [2000]. They 
have calculated the trajectories of the pickup ions in the 
prescribed magnetic field obtained from their earlier MHD 
simulation. 

There have been two most recent simulations one in MHD 
and the other a hybrid simulation. The first is a multi- 
species MHD simulation by Nagy et al [2001], which 
separately accounts for three generic species, light (e.g^., H , 
H 2 \ H 3 h ), medium (e.g., N + and CH 3 ) and heavy (N 2 \ and 
HCNH + ). This permits the effects of exospheric mass 
loading, chemisca! reactions and ion-neutral collisions to be 
treated. It reproduced the global features (magnetic barrier, 
magnetotail and distributions of the major ionospheric 
species) and implied an escape flux of major ionospheric 
(heavy) species to be 6.5 X 10 24 s’ 1 . The other is a full hybrid 
simulation, like those discussed by Brecht [chapter II.3], 
which has been performed by Brecht et al. [2000]. The 
hybrid simulation accounts for particle ions and fluid 
electrons in addition to computing the magnetic field. It 
includes a similar grouping of species classes. Although this 
approach is better justified for Titan, it is also unfortunately 
more computationally expensive, so as with the case of lo 
above, there are always trade-offs in comparing one 
simulation method to another. 

The hybrid simulations of Brecht et al. [2000] quite 
naturally produced an asymmetic plasma flow around Titan, 


something that is not possible with the ideal MHD 
simulations. This is produced by the inclusion of Hall 
terms enabled by the explicit inclusion of finite gyroradius 
microphysics. However, the general scale of the larger 
disturbance created by Titan’s ionosphere and mass-loading is 
actually quite similar to the MHD simulations. 


4. FINAL COMMENTS 

The title of the paper by Cloutier et al. [1978], which we 
have cited elsewhere in this review, is "A cometary 
ionosphere model for Io." This clearly portends the 
subsequent scientific and methodological connections 
between planetary satellites and comets in the space science 
literature. In fact the reference sections of papers in both 
fields show clearly that each field benefits from the 
contributions of observers, experimentalists, and theoreticians 
from the other. Here we have presented a review of the many 
common properties, features and modeling approaches of the 
tenuous atmospheres, ionospheres, and outer radiation, fields 
and particles environs of comets and outer planet satellites. 
Although much progress has been made in understanding 
global structures and local measurements, there is much work 
to do in terms of fundamental theoretical modeling of the 
combination of global morphologies with local microphys- 
ics. As computation power continues to advance during the 
near future it will be more and more possible to incorporate 
complicated detailed physical processes in our models. The 
challenge will be to stay at the forefront, constantly 
stretching computational capabilities, in order to improve our 
understanding of these fascinating subjects. 
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